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Based on the results of a previous paper (Paper I), by performing the geometrical mapping 
via coherent states, phase transitions are investigated and compared within two algebraic cluster 
models. The difference between the Semimicroscopic Algebraic Cluster Model (SACM) and the 
Phenomenological Algebraic Cluster Model (PACM) is that the former strictly observes the Pauli 
exclusion principle between the nucleons of the individual clusters, while the latter ignores it. From 
the technical point of view the SACM is more involved mathematically, while the formalism of 
the PACM is closer to that of other algebraic models with different physical content. First- and 
second-order phase transitions are identified in both models, while in the SACM a critical line also 
appears. Analytical results are complemented with numerical studies on a-cluster states of the 20 Ne 
and 24 Mg nuclei. 
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I. INTRODUCTION 

In a former contribution called from now on Pa- 
per I, we investigated the geometric mapping of the Phe- 
nomenological Algebraic Cluster Model (PACM) and the 
Semimicroscopic Algebraic Cluster Model (SACM) @,[|[. 
The first does not observe the Pauli exclusion principle 
while the second one does. The PACM belongs to the 
same family as models like the vibron model In 
both type of models we considered the same Hamilto- 
nian, while the model space of each is quite different. 
In the SACM, observation of the Wildermuth condition 
Q means that the number of relative oscillation quanta 
is restricted from below while in the PACM the relative 
oscillation quanta start from zero. 

The method of geometrical mapping is gleaned from 
II for the SACM, which reduces to the usual one @-[|| 
when no Pauli exclusion principle is taken into account. 
We showed in [l[ that the differences in the mapped geo- 
metrical potential are large. However, within the PACM 
one can reproduce the geometric potential of the SACM 
by including very complicated higher-order interactions. 

We also included in [l[ a discussion on how to define 
correctly the coherent state parameter, such that for a 
large total number of bosons (N, which plays the role 
of a cutoff) the mapped potential will be independent of 
the cutoff N. 

In this paper we concentrate on the study of phase 
transitions and show that not only second-order phase 
transitions may occur, but also first-order phase transi- 
tions. Furthermore, in the SACM a critical line appears 
beyond which no phase transition occurs. Thus the struc- 
ture of the phase diagram of the SACM will be much 
richer. 

We will also apply the mapping to two kind of systems, 



one with two spherical clusters and the other one with 
a deformed and a spherical cluster. Numerical studies 
are added. We will show that the PACM leads to incon- 
sistencies when applied to real nuclei, while the SACM 
performs well. 

The paper is structured as follows: in section |H] a gen- 
eral discussion on phase transition, both in the PACM 
and SACM, is given, independent of a particular system. 
In Section [ill! results are illustrated with two particular 
cluster systems, the a-cluster states of 20 Ne and 24 Mg. 
In the former, both clusters are spherical, while in the 
latter, one of them is deformed, leading to a more com- 
plex physical situation. Finally, in Section llVl conclusions 
are drawn and a discussion is presented on the differences 
between the PACM and the SACM and their importance 
in the study of nuclear clusters. 



II. PHASE TRANSITIONS IN ALGEBRAIC 
CLUSTER MODELS 

In what follows we apply the formalism developed in 
Paper I [TJ] to discuss phase transitions in the SACM and 
PACM. Before that we present the general framework 
within which the discussion will be implemented. 



A. Definition of a phase transition 

Phase transitions are investigated u sing the following 
steps and the recommendations of [10l - tl2l ]. This method 
can be applied to any system and does not depend on the 
notion of symmetries. This presentation does not need 
the language of catastrophe theory [l3[, which simplifies 
considerations. 
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a) In the first step the minima of the potential energy 
surface (PES) are determined in the space of the collec- 
tive variables a rn . In the present case there is only one 
relevant variable a. This is due to the fact that the dis- 
tance vector between the clusters can always be aligned 
along the intrinsic z-axis which connects the two clusters. 
The extrema are obtained from 



dV_ 

da 



= 



(1) 



This determines the position of the extrema at 
oti (i = 1,2,..) and the values V(c<i) of the potential 
there. The oti are the values of the variable a at the i th 
minimum, which is a function of the interaction param- 
eters p = (pk), with pk as a short-hand notation for the 
parameters k = 1, 2, rife, with being the number of 
parameters. 

b) Once the extrema are obtained, one determines at each 
minimum the first and second derivatives of the potential 
with respect to the parameters of the model, i.e., 



d n V(ai 



1,2 



(2) 



When we discuss the PACM and SACM further below, 
we will give first a general discussion on the properties of 
phase transitions, independent of the values of a partic- 
ular subset of interaction parameters. In the subsequent 
concrete applications, however, we will fix most interac- 
tion parameters and vary only x and y, which control the 
transition from one effective symmetry (l4| to another. 
This will give us particular curves in the space of phase 
transitions. 

c) The locations of phase transitions are determined by 
identifying the borders in the parameter space where at 
least two minima are at equal energy, i.e., 



V(c 



(3) 



for some index values i\ and «2 of the minima. This 
results in a relation 



f(pi,...,Pn p ) = 



(4) 



between the parameters, which allows one, in principle, 
to express one parameter in terms of the others (see the 
particular examples further below). 

d) The phase transition is of order m when, up to n = 
m — 1, the derivatives of the potential with respect to the 
free parameters are equal at the point of phase transition, 
while the m th derivative of the potential, with respect to 
its parameters, is discontinuous at the point of the phase 
transition, i.e., 



d n V(a n ) d n V(a. l2 ) 



dpi 
dp^ 



dpi 

d m V(a i2 ) 



for n < m 



(5) 



Note again that this procedure for determining the or- 
der of a phase transition is quite general and does not 
depend on identifying the phase with a dynamical sym- 
metry. Thus, this procedure can also be applied to sys- 
tems that do not exhibit dynamical symmetries. 



B. A possible explanation for the differences in 
using coherent states and a numerical approach 



In the numerical calculations further below, we will 
see that, though the use of coherent states shows a clear 
phase transition, this is not the case in the numerical 
study, because the total number of bosons N is fixed. In 
what follows, we indicate why this is so. 

The main reason for using coherent states is that the 
phase transition is given by the crossing of two poten- 
tial minima. In the case of a first-order phase transition, 
the two minima are well separated by a barrier and a 
simple crossing takes place. For the case of a second- 
order phase transition, the spherical minimum coincides 
with the "deformed" one at a = at a certain point in 
the parameter space. Once one minimum is deeper than 
the other one, an immediate jump by hand to the global 
minimum is made, using coherent states. This appears 
to produce the well defined phase transition. However, in 
numerical calculations within a finite system, the states 
are a mixture of states in both minima. Around the point 
of phase transition, the eigenstates of the Hamiltonian 
are a mixture of states in both minima. Only far away 
the states are located clearly in one or the other mini- 
mum. This produces, for finite N, a gradual change in 
the wave function and as a consequence a gradual change 
of the control parameters, like the energy or the expecta- 
tion value of n„, which define the phase transition. For 
N — > oo, the discontinuities in the derivatives of the po- 
tential will be clear cut. 

As a consequence, the apparently clear phase transi- 
tion for any number of N is not a real one, because N is 
finite. Nevertheless, the use of coherent states shows the 
points of phase transitions and of which order they will 
be for N — > oo. 



C. Study of phase transitions in the SACM 

Paper I [l[ contains the general expression of the 
Hamiltonian (in subsection 2.2) and the geometrically 
mapped SACM potential (in subsection 4.1). Here we 
recall only the essential formula, Eq. (30) of Paper I, 
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necessary for the present discussion: 



(H) = C(x,y)-(b + b)xy[A(x,y)a 2 ^^ 



F 00 (a 2 ) 



-B{x,y)c 



-C(x,y)a' 



,6^33 (a 2 ) 



^22 (a 2 ) 
F 00 (a 2 ) ' ~ F 00 (a*) 

; <- 2 (« 2 ) \ 

^00 (« 2 ) J 



(6) 



We also note that the x = case relevant to the 5*0(4) 
to 50(3) phase transition has to be discussed separately. 
(See Eqs. (35) to (38) in Paper I Q.) 



1. Phase transition diagram for the SACM: General 
considerations 

The complex structure of the geometrically mapped 
potential complicates an analytic treatment of the prob- 
lem. It is therefore essential to formulate a set of criteria 
facilitating a straightforward way to determine the order 
of phase transitions. 

Equation demonstrates that the explicit depen- 
dence of the potential on the parameters A, B and C 
is linear of the type 



V 



E 



p k a m »f k (a) 



(7) 



with m,t > 1 and pk being a short-hand notation for the 
k th parameter. The fk are given by ratios of the functions 
Fpq{a) and are always greater than zero for a ^ 0. Only 
the function FJ^~ 2 approaches zero for a — > oo. (Sec 
Eq. (39) in Paper I QJ.) 

According to the general discussion on phase transi- 
tions in Subsection III Ai let us now turn to the potential 
minima located at on. We investigate their dependence 
on the parameters pk, standing for A, B and C. The 
structure of (JTJ) guarantees that ct\ = is always an 
extremum. Furthermore, not only its first-order deriva- 
tives, but the function value is also zero at ct\ = 0. It 
is thus sufficient to focus on the second, deformed min- 
imum, with «2, for which the following consideration 
holds. There are two possibilities: 

a) a 2 > 0, or 

b) a 2 = 0. 

The consequences are seen by determining the first ab- 
solute derivative of the potential with respect to the pa- 
rameter pk- This first derivative is given by 



dV 
dpk 



dV + dV da t 
dp k dan dp k 
dV 
dpk 



(8) 



because j^- vanishes at the minimum. Taking into ac- 
count (0, this is further expressed as 



dV 
dpk 



- — = a 2 nk f k {a 2 ) 



(9) 



For case a) this expression is different from zero (re- 
member that 6t2 > 0), while the derivative within the 
spherical minimum is equal to zero. In other words, there 
is a first- order phase transition. 

For case b), © is equal to zero, and therefore the phase 
transition must be of higher order. To determine which 
order, second-order derivatives arc also needed: 



d 2 V 
dpj 




da™* 
dpk 



fk{a 2 ) 



dpk 



a 2 



(10) 



Here we used the fact that in © there is no explicit 
dependence in pk left. Since in this case a.2 = 0, (fT0|) 
reduces to 



d 2 y 
dpi 



da2 k 
dpk 



fk(a 2 ) 



(11) 



As both the fk functions and the partial derivatives are 
different from zero, so is (|TT|) in general, i.e. the phase 
transition is of second order in case b). 

In conclusion we have the simple identification of phase 
transitions: 

a) When the deformed solution 0,2 is different from zero, 
then there is a first-order phase transition. 

b) When the deformed solution 0,2 is equal zero, then 
there is a second-order phase transition. 



2. General discussion of phase transitions in the SACM 

In Fig. [1] the surface where a phase transition takes 
place is plotted in the space of the independent param- 
eters A, B and C. The solid line marks a change from 
first-order to second-order transitions. For larger, posi- 
tive C the transition is of second order, while for smaller, 
negative C it is of first order. For completeness we show 
in Fig.[5]the phase space diagram for the case when x = 0, 
as it is the case for the 5*0(4) to 50(3) transition. 

A remarkable finding is that at approximately 
C ~ -15 the surface of phase transition ceases to exist. 
Beyond that point, no phase transition can be observed, 
i.e., the straight dashed line in Fig. Q] at approximately 
C ~ — 15 represents a critical line. This is shown by a 
dashed line in Fig. [T] One then can trace a straight line 
from below the surface, around the critical line, ending 
up above the surface without passing through a phase 
transition, which is similar to the critical point in the 
two-dimensional phase diagram of water. Fixing all in- 
teraction parameters except x and y (as we will do in 
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the numerical applications) , passing from one dynamical 
symmetry limit to another one will trace a line in this 
three-dimensional space. Depending on the fixed param- 
eters, this line will or will not cross the surface of phase 
transition. 





FIG. 3: The variable 612 of the deformed solution, as a func- 
tion in A and C. B is fixed by the requirement that one is 
at a point of a phase transition. The solid line marks the 
transition from a second- to a first-order phase transition, the 
dashed a 'critical line'. 



FIG. 1: The phase space diagram of the SACM as a func- 
tion of the independent parameters A, B and C. The solid 
line marks the change from a second- to a first-order phase 
transition, while the dashed line represents a 'critical line'. 



- Phase Diagram for the case when x=0 for every point 




-0.15 -0.10 

IT 



FIG. 2: The phase space diagram as a function of the inde- 
pendent parameters A and B for the SO(4) to SO(3) tran- 
sition, within the SACM. As defined in Paper I, within the 
SO (4) dynamical limit the x = and one has to define a new 
representation of the potential, involving the new parameters 
A and B. (See also Paper I [l[ for more.) 

Figure [3] displays 012 corresponding to the deformed 
solution at the point of phase transition as a function in 
A and C . (Remember that the spherical solution corre- 
sponds to the always existing extremum at a = in the 
case it is a local minimum.) Each point of the surface 



represents also a given B at which the phase transition 
occurs, i.e., B is fixed by the requirement that there is 
a phase transition. In this figure the solid line also rep- 
resents the change from one type of phase transition to 
the other one. For larger, positive C 612 = holds, and 
according to the discussion of subsection III C 1[ it corre- 
sponds to a second-order phase transition. For smaller, 
negative C the 0.2 > holds at the point of phase transi- 
tion, i.e., it corresponds to a first-order phase transition. 

The whole phase structure of the system is illustrated 
in Fig. [U which corresponds to the fixed value A = 10. 
The solid line represents a cut through the phase transi- 
tion surface. For C > the phase transition is of second 
order, while for C < is of first order. The Roman 
numbers indicate the following regions: i) Region I cor- 
responds to the existence of two minima (one spherical 
and the other deformed), with the deformed one as the 
global minimum, ii) Region II also corresponds to two 
minima with the spherical as the global one. iii) In region 
III there is only one spherical minimum while iv) in re- 
gion IV there is only one deformed minimum. The region 
denoted by a zero refers to potential with no minimum. 
The other dashed lines do not indicate phase transitions, 
rather they separate the areas where two minima exist 
and the ones where only one minimum exists. The line 
of phase transition ends at approximately C f=s —15, cor- 
responding in this case to a point. 

There is another point on the solid line, where the 
first-order phase transition turns over into a second-order 
phase transition. A distinction has to be made between 
this and the critical point of the phase diagram of water, 
where before the critical point the phase transition is 
of first order, while at the critical point it is of second 
order. In our case, however, the phase transition does 
not end at this point, rather it continues as second-order 
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one in the domain with C > 0. One can now draw a 
curve which encircles this line, going from a spherical 
phase to a deformed one. The value of C at which this 
phenomenon happens, varies with A and B. 




A=10,C = -10 




FIG. 6: (color online). The value of ai at A = 10 and 
C = — 10 as B moves through regions. 



FIG. 4: (color online). A cut through the phase at A — 10. 
The notation is explained in the text. 




FIG. 5: (color online). Several potentials for different values 
of B and C, for A = 10 are shown. Some potentials approach 
a positive value, i.e., for N — > oo they goes to + oo. Others 
approach a negative value, i.e., for N — » oo they approach 



Fig. [5] gives a sample of energy functionals at a fixed 
N (tiq = 8 and N + no = 20) in different regions, the 
specific points being marked with lower-case letters in 
Fig- HI The upper left panel shows functionals either side 
of the first-order phase transitions, points a) and b). The 
upper right one shows graphs either side of the second- 
order phase transitions, points c) and d). The lower left 
panel shows a functional either side of the dotted line in 
the Region III of greater B, points e) and f), showing 
that while each have a barrier, one is bound as a — > oo, 
and the other is unbound. These tails arc dependent 
on N, with the unbound approaching infinity and the 
bound negative infinity as N — > oo. Thus, this region 
is unphysical. The Region III of lesser B, shown in the 
second panel, has no barrier, and is bound as a — > oo. 
The fourth panel shows a fully unbound potential with 
no minima from Region 0, point g). 

For completeness, Fig. |H] traces the value of ct2 for 
14 < B < 24 at C = -10 in Fig. 0J which thus represents 
a movement between several regions in that figure (and is 
therefore analogous to being a subset of Fig. [3]). It shows 
that the critical point corresponds to ct2 asymptotically 
approaching oo at B slightly greater than 22. 

Note that the discussion of phase transitions is com- 
pletely independent of the system considered, which is of 
great advantage. In the next section we will consider par- 
ticular systems and study the properties of their phase 
transitions. 



The horizontal dotted line represents a division (see 
discussion in III A[) : above this line the potential ap- 
proaches a negative infinite value for N — > oo, i.e., the 
potential gets unstable and the cluster system dissolves, 
once getting to the "deformed" solution. Below that line 
the potential approaches plus infinity for — > oo, thus 
the deformed solution corresponds to a stable cluster sys- 
tem. 



D. The PACM case: Pauli principle not taken into 
account 

Here we apply the results obtained in subsection 4.2 
of Paper I [lj. We display only the expression of the 
normalized potential in Eq. (46) of Paper I: 



V 



{A/3 2 - Bf3 A + /3 e } 



(12) 
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We also remind the reader that the SO(A) to 50(3) tran- 
sition needs a separate discussion, because in that case 
the potential reduces to a quartic one. (See the discussion 
in subsection 4.2 of Paper I [l[.) 



1. Phase transitions 

The extrema are found by setting the first derivative 
with respect to (3 equal to zero. We obtain, in general, 
two solutions 



solution 1 : /3i = 



(13) 



noting that in this case the potential value is always 
V((3i) = 0, and 



solution 



2: Pi = \(b±^B 2 ~ZA) , (14) 



where the value of the potential at the extremum with 
the positive sign is given by 



6A- B 2 - B\J B 2 - 3A) . (15) 



The barrier maximum is obtained by choosing the nega- 
tive sign in front of the square root expressions. The last 
equation is acceptable only when the solution in ([141) is 
real and positive. The solution is real for A < B 2 /3. 
When B < and A > holds, no deformed solution ex- 
ist. In fact, inspecting Eq. (fT2j) for A > there is always 
a minimum at ,9 = 0, while a negative B always gives a 
positive contribution, +/3 4 . When A and B are negative, 
the potential is negative at small /? and is turned over by 
the positive contributions ~ +/3 4 and +(3 6 . In this case, 
only a deformed minimum exists. 

The structure explained above is summarized in Fig. [7] 
The horizontal line is the B axis while the vertical one 
is the A axis. On the left hand side (B < 0) a spherical 
minimum exists for A > 0, while for A < the minimum 
is deformed. The phase transition takes place at the line 
A = 0. This will be discussed further below as case b). 
For B > and A < 0, a deformed minimum always 
exists. For A > and below the dashed curve A = B 2 /3 
a deformed minimum coexists with a spherical minimum 
until A = is reached. The solid curve is determined 
by requiring that the potential minima at /3i and /?2 are 
degenerate, i.e., 



which leads to the condition 



(16) 



(17) 



Between the upper dashed and the solid curves the spher- 
ical minimum is the global one, while below the solid 




FIG. 7: (color online). The parameter phase diagram for the 
PACM [l|, m. The horizontal axis is B, while the vertical 
axis corresponds to A. In Region I and Region II two minima 
exist, one spherical and one deformed. In Region I the global 
minimum is the deformed one, while in Region II it is the 
spherical minimum. In Region III only a spherical minimum 
exists and in Region IV the only minimum is a deformed one. 



curve the global minimum is the deformed one. Crossing 
the solid line a phase transition takes place, which will be 
called below case a). It is of no surprise that the phase 
diagram in Fig. [7] is similar to Fig. 2 of (l7j. 

Case a), B > 0: 

The potential is exactly zero at /?i =0, and thus 



<9"V(/3i = 0) 
dA™ 



(18) 



for all ri > 0. Contrary to this, the first derivative of 
the potential at /?2, the deformed minimum, taken at 
A = B 2 /4 where the height of the deformed minimum 
equals that of the spherical one, is given by 



DA 



l -B 

2 



(19) 



which is clearly different from zero for B > 0. Thus the 
phase transition is of first order for B > 0. 

We can repeat the calculation for the derivatives in B, 
leading to the same conclusion. Also, any combination 
of curves given by the variable (cA — sB) (with c = cos0 
and s = sin0, with <fi as the angle between the tilted 
straight line and the A axis) leads to the same conclusion. 

Case b): B < 

We will vary A from positive to negative values for a 
fixed B. The same steps as taken for case a) can be ap- 
plied, setting B = —\B\. The first and second derivatives 
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of the potential with respect to A are given by 

dA 3 V V VB 2 - 3A 

1 



<9,4 2 



2VS 2 - 3A 



Evaluating these derivatives at A = and B 
gives 



9.4 



= 



d 2 V 
dA 2 



1 

25 



(20) 



\B\ 



(21) 



For this case we also have to determine the derivatives 
with respect to B, given at the point of phase transition, 
i.e., setting A = afterward, and B = — \ B |, 



parameter was not fixed, but rather it was chosen accord- 
ing to the harmonic oscillator constant corresponding to 
the unified nucleus. The fixed interaction parameters are 
determined in such a way that in the dynamical symme- 
try limits the spectrum appears with the same scale as 
the physical measured one. The spectrum of the real nu- 
cleus would probably correspond to a single point in the 
(x, y) parameter space. However, our aim is not repro- 
ducing the exact spectrum, rather to reach a conceptual 
understanding of phase transitions when going from one 
dynamical symmetry to another, and investigating the 
conjecture that a phase is defined by an effective sym- 
metry. One alternative method would be to adjust sev- 
eral cluster systems and to try to find a series of systems 
which, for example, would cross the surface of phase tran- 
sitions at one point. We do not follow this way, rather 
postponing it for later consideration. 



dV 


1 - 


dB 


~9 . 


d 2 V 


2 ■ 


dB 2 


~9 . 



-\B\ + ^/\B\ 2 -3A 



A. Two spherical clusters: ie O+a 



J Ne 



\B\ 







V\B\ 2 -3A 
(22 



Changing B corresponds only to moving on the line of 
phase transition, which does not make sense, because one 
has to cross the line of phase transition. However, a com- 
bination of changing A and B corresponds to an inclined 
straight line, crossing the line of phase transition. This 
means that one has to consider a general line given by 
the variable C = cA — sB, with arbitrary c = cosc/> and 
s = sine/), with cj) being an arbitrary rotation angle. De- 
termining at the point of phase transition the first- and 
second-order derivatives of V with respect to C leads to 



dV_ 

dC 
d 2 V 
dC 2 



-> c 



dV 
dA 
,d 2 V 



-> c 



dA 2 



(23) 



Because the first derivative of the potential with respect 
to A is zero, while the second is different from zero at 
the point of phase transition, we can claim also that the 
first derivative with respect to C is zero, while the sec- 
ond is different from zero. We conclude that the phase 
transition is of second order. 



III. NUMERICAL STUDIES 

This section deals with two widely known cluster sys- 
tems: 20 Ne as 16 0+a, where both clusters are spherical, 
and 24 Mg as 20 Ne + a, where one of them is deformed. 

Since our aim is studying transitions from one particu- 
lar dynamical symmetry limit to another one, we first fix 
all the interaction parameters except x and y. The fun 



In this case the only degree of freedom is the radial 
motion, as the clusters do not have an internal structure 
apart from the fact that they are composed of fermions, 
which have to obey the Pauli exclusion principle. There- 
fore the cluster representation is (AcA'c) = (0,0). As 
described in Paper I the 50(3) limit does not exist 
as an independent limit in this case (the 50(3) Hamilto- 
nian is a reduced version of the the SU (3) Hamiltonian) , 
so the y = 1 choice has to be made. The only transition 
to consider is thus between the 5t/(3) and the 50(4) 
limits. 

Concerning the determination of the parameters, one 
has to take into account that some parameters appear 
in both dynamical symmetry limits, like 7 preceding L 2 
(see Eqs. (10) and (11) in Paper I [l|). We first deter- 
mine the parameters in the SU(3) dynamical symmetry 
limit, which fixes 7, and then we determine the remain- 
ing parameter c, which appears in the 50(4) dynamical 
symmetry limit. The terms C2(Ac, fic) an d L c do not 
contribute to the Hamiltonian in this case, so the cor- 
responding parameters are kept zero. Note that in this 
case L 2 R = L 2 . 



1. The SACM 

In the first step the parameters are adjusted within 
the SU(3) limit, setting x = 1 and y = 1 and in the 
50(4) limit, setting x = and y = 1. The parameters 
are depicted in Table HI The SACM yields reasonable re- 
sults, because the ground-state band belongs to n n = 8, 
(A, fj,) = (8, 0), where n T = 8 corresponds to the minimal 
number of relative oscillation quanta no required by the 
Wildermuth condition. The first excited + state cor- 
responds to a 2hiv excitation and naturally lies at high 
energy as required by the experimental data. The spectra 
in the 50(4) and SU (3) limits are shown in the left and 
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TABLE I: Parameter values defining the 16 0+a interaction. 
See Eq. (11) in Paper I [l|. 

Hamiltonian 



a a 

-0.500 0.000 

c ac 

0.250 0.000 



7 

0.208 

„(!) 



dClus 

0.000 

t 

0.000 0.000 



b b 

0.000 -0.009 



Clusters 


Al fj,! 


iVo.i 


ft A 2 


/U2 


Ao,2 


ft 





0.00 


0.00 





0.00 





Quanta 






no 




N 




13.2 




8 




12 





right extreme of right panel of Fig. [TT} respectively. The 
spectrum of experimental 20 Ne states, each correspond- 
ing to this clusterisation, is shown in the left panel. As 
already mentioned, the real nucleus will lie somewhere 
between x = 1 and x — 0, though, the fit at x = 1 is 
acceptable. 




x can be seen in Fig. |H1 In summary, the situation in 
this example corresponds to a phase transition of second 
order. 




SU(3) 



FIG. 9: (color online). The parameters A, B and C as a 
function of x, considering the 5(7(3) to SO(4) transition. The 
inset shows a region of large change in A, B and C with change 
in i. A selected range of values is shown in the inset, for a 
better reading. 




SU(3) 



FIG. 8: (color online). Part of the SACM phase space dia- 
gram, where the curve of the parameter values, as a function 
of x, crosses the surface of phase transition. The crossing 
occurs at approximately x = 0.6 and happens at the surface 
related to the second-order phase transition. 



FIG. 10: (color online). The expectation value of (vertical 
axis) as the function of x. Starting from x — 1 the expectation 
value is 8, given by the lowest possible number of n quanta. 
From x — 0.6 on we observe a rise in the expectation value, 
reaching n v = 11 at x — 0. The phase transition takes place 
at x — 0.6. 



In order to see if a phase transition appears, and of 
which order it is, we added a curve to the (A, B, C) phase 
space in Fig. [5J depicting the transition from x = 1 to 
x = 0. The figure shows only the relevant part of the 
phase space, i.e. the one where the curve crosses the sur- 
face associated to the second-order phase transition. This 
occurs at approximately the x = 0.6 parameter value. 
The actual values of A, B and C associated with this 



Figure [TO] displays how the expectation value of 
changes as the function of x. The lighter (yellow on- 
line) and darker (red online) curves depict the result of 
the geometrical mapping and the numerical diagonaliza- 
tion, respectively. From x = 1 up to the point of phase 
transition, the effective SU(3) limit is realized and the 
expectation value is equal to the minimal number of tt 
bosons, i.e., no = 8. Below x = 0.6 the darker (red) 
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FIG. 11: (color online). The lowest energy states as the 
function of x for 16 0+a — > 20 Ne. The spins and parities are 
depicted in the legend. 

curve begins to rise slightly, indicating that the structure 
of the system is changing. 

The right panel of Figure Qj] shows the lowest energy 
states of this clusterisation, those fitted, as the function 
of x. A smooth behavior is observed, with an accumula- 
tion of states near x = 0.6. This is in accordance with 
the expectation value of in Figure [TO] 

2. The PACM 

This type of model was considered in [l8j , where a pure 
schematic investigation on possible phase transitions was 
presented. The model space was restricted to low ex- 
citations. No parameter fit was applied to a physical 
system. We show here that the relation to a physical 
system is of utmost importance and can discriminate be- 
tween physical and non-physical models. 

As a first step, we tried to adjust the parameters of 
the model in the SU(3) limit (x = 1 and y — 1) and 
in the SO(4) limit (x = and y = 1). However, we 
already encountered severe problems in the SU(3) limit: 
The model space for even angular momentum starts with 
n n = 0, thus the lowest states are comprised by ra, = 0, 
(A, /i) = (0,0), which contains one J 71 " = + state, and 
by «7r = 2, (A, [i) = (2, 0), which contains a + and a 2+ 
state. The problems encountered are: 

i) Since the factor in front of L 2 has to be positive 
and supposing that the ground state belongs to n T = 
and the next excited positive parity state belongs to 
n w = 2, the first excited + state will always be lower 
than the first excited 2 + state. This contradicts the 
experimental spectrum with E(2f) = 1.634 MeV and 
E(0+) = 8.7 MeV. 

ii) Since the 2f state belongs to a 2hoj excitation, with 
huj = 13.2 MeV, the quadrupolc-quadrupole interaction 



has to be unnaturally strong in order to shift the energy 
to 1.634 MeV. This, in turn, will move very high ex- 
citations to low energy, even below the supposed ground 
state with = 0. 

iii) Due to the completely different SU(3) structure of 
the states within the ground state "band", one can not 
talk about a rotation band. 

Restricting to the SU(3) limit, the eigenvalues of the 
Hamiltonian are given by 

E = hum n + (a - bn^rinirin + 3) + jL(L + 1) . (24) 

Here we already see that for 7 > and a fixed n T , higher 
spin states are higher in energy. In order to adjust the 
2^ and the 0^ states to the experimental energies, the 
fitting routine assigns to both the ground state and the 
2+ state a different n^. Using N = 20, as a result we 
obtain = 20 (the total number of bosons was set to 
be 20) for the 0^~ and the 2+ state, while the 0^ state 
belongs to = 0. When we change the total number 
of bosons, we get similar results. The mere fact that we 
have to involve states with = N indicates that no 
convergence is achieved, considering that N represents a 
cut-off value. 

Similar results are also obtained when the SO (4) limit 
is considered. In the 50(4) limit, the energy is given by 

E = ^{N -uj){N + uj + 2)+ 1 L{L + 1) , (25) 

where uj refers now to the 50(4) quantum number with uj 
= N, N — 2, or 1. The lowest state is normally taken 
as uj = TV, which contains L = 0, 1, ...,lu. Thus, choosing 
N = 20, the ground-state band (u> = 20) is composed of 
the angular momentum states L = 0, 1, 20 and the 
first excited band with even spin (uj = 18) is given by the 
states L = 0, 1, 18. Thus, the first excited + state 
can be set at higher energies than the first excited 2 + 
state, adjusting the parameter c. Everything seems to 
be in order, except for the problems which the following 
discussion demonstrates. 

The difference with respect to the SU(3) limit is that 
a state is a mixture of many basis states in SU(3). We 
adjusted the spectrum of 20 Ne to the 50(4) limit and 
confirmed that for L = the expectation value of the 
operator n, is given by [l|| 

, > N-l 

K) = ■ (26) 

The problem here is that when the cut-off is increased, 
the structure of the states change: (n n ) increases, im- 
plying that no convergence has been reached. This also 
implies high shell excitations, if we assume that the two 
clusters are moving in a shell model mean field, which 
has been proven in many microscopic calculations j20| . 
In the PACM, however, the mean field hujn^ has no spe- 
cific meaning, i.e., the parameter Tiuj — > d can be very 
small. As shown above, the spectrum can be easily ad- 
justed within the £0(4) limit. The fact that it adjusts 
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TABLE II: Parameter values for the 20 Ne+a interaction. See 
Eq. (11) in Paper I 0. 

Hamiltonian 



a a 7 

-1.396 -0.136 0.197 



c ac 
0.470 0.079 



dClus 

0.000 

t 



b b 
-0.116 0.045 



0.053 0.664 



Clusters 



Ai 



Mi N o,i Pi 
48.5 0.73 



A 2 




//2 





A*0,2 

4.5 



ft 




Quanta 



12.6 



8 



JV 
12 



the spectrum suggests that there must be some truth in 
it. Nevertheless, the basic degrees of freedom (clusters 
plus relative motion) cannot be interpreted as real clus- 
ters or relative motion, as is done in microscopic cluster 
studies, but must be in a complicated relation with them. 
Just what relations these are remains a big problem. 

This result demonstrates that the model shows incon- 
sistencies and the reasons are exposed in the arguments 
i) and ii) above. The intcrprctational problems of the 
PACM have already been indicated in jl4j |. 

A possible solution to this problem is to redefine the 



pairing operator as 



rtl 



R 2 a 



i.e., the in- 



troduction of a new parameter R 2 which can be set pro- 
portional to 1/N. In this way, the N dependence can 
be eliminated and physical results can be expected. This 
procedure was adopted in (2lj . where the vibron model 
was extended to three clusters describing the 12 C nucleus 
as an oblate symmetric top. Though, in [2l| the SO (7) 
limit in a U(7) algebraic model, which describes three 
clusters, is considered, the same ideas can be translated 
here. 




i i i i i i i i i i i i i i i i i i i i i i i i 

0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.8 0.6 0.4 0.2 

SO(3) y SO(4) x SU(3) y SO(3) 



FIG. 12: (color online). Expectation value of for the 
three transitions i) SU(3) to 50(4), ii) 5*17(3) to 50(3) and 
iii) 50(4) to SO{3). 




exp- 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.8 0.6 0.4 0.2 



SO(3) 



SO(4) 



SU(3) 



SO(3) 



B. One spherical and one deformed cluster: 

20 Ne+a 24 Mg 

This is the first example where the cluster part has a 
structure owing to the deformed 20 Ne. 



1. The SACM 

In order to analyze the transitions between the SU(3), 
50(4) and SO(3) limits, first the parameters were fixed 
in the three limits. The results are displayed in Table [TT1 

Figures Q2] and Q2] display the expectation value of n n 
and the lowest states in energy, respectively, for all three 
transitions. The darker (red online) curve in Fig. 1121 
shows the results of the numerical diagonalization, while 
the lighter (yellow online) curve shows the one of the 



FIG. 13: (color online). The lowest states in the three tran- 
sition i) 577(3) to 50(4), ii) 5(7(3) to SO (3) and iii) SO (4) 
to SO{3). 



geometrical mapping. The results are qualitatively sim- 
ilar to those encountered in the former example. The 
expectation value of n„ starts in the SU(3) limit at 8 
and increases towards the SO(4) limit. The transition 
is smooth in the numerical calculation, while in the geo- 
metrical mapping the transition is pronounced well. The 
transition is indicated by a sudden change in the slope 
at above x = 0.2 in the central panel of Fig. [T3J The 
energy spectrum in Figure [T3] does not show a particular 
structure at points of phase transition. Therein, the left- 
most panel is again the experimental spectrum used, and 
the right three are theoretical results as x and y are ad- 
justed. In the SU (3) to 5*0(3) transition no phase change 
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appears, because the global minimum of the potential is 
always at a=0. One can observe a distinct behavior as 
a function of x and y between the states with positive 
and negative parity, marked as solid and dashed curves, 
respectively. The latter are more sensitive to the change 
in x and y. 

Similar properties can be observed in the other limits, 
i.e., see the leftmost panel in Fig. Q21 In the rightmost 
panel of Fig. [T2] no phase transition occurs. The energy 
spectrum also shows an accumulation of states at low 
energy at the point of phase transition. 

2. The PACM 

Here we find similar inconsistencies with respect to the 
model space as in the case of 20 Ne in IIII A 2l The cluster 
irrep of 20 Ne is (8,0), while the relative oscillation irreps 
are ( 71^,0), with n w = 0,1,2,... Restricting to positive- 
parity states only, the lowest energy model space, in the 
SU(3) limit, consists of (8, 0) at Ohui, and (10,0) and 
(6,2) for the 2Huj excitation. The ground-state band is 
a (8,0) irrep and the lowest K = 2 band is the (10,0) 
irrep at 2Tku. Even if we change Tiuj to an arbitrary small 
parameter a\, in order to bring the 2fiui states down in 
energy, the internal structure of the ground state band 
is not what we expected, namely (8,4). Also there is no 
B(E2) transition between the ground state band and the 
K = 2 band, because they belong to different irreps. (See 
the discussion in Subsection 2.1 in Paper I [J.) 

For the 50(4) limit one obtains a satisfactory fit, but 
again with the fact that the expectation value of the num- 
ber operator n, depends on the cut-off N. 

For these reasons, we do not present figures of the spec- 
tra and expectation values. 

IV. CONCLUSIONS 

Phase transitions were investigated in two algebraic 
cluster models, one of which observed the Pauli exclusion 
principle between the nucleons of the individual clusters 
(SACM), while the other (PACM) did not. This analysis 
was based on the results of a previous work [l[ , in which 
the geometric mapping of the two models had been per- 



formed using the coherent state formalism, leading to ap- 
propriate potential energy surfaces. The dynamical sym- 
metries of the SACM and PACM had also been identified 
in [l[, and in the present analysis special attention was 
paid to transitions between phases associated with the 
SU(3), 50(3) and 50(4) dynamical symmetries. The 
potential energy surfaces depended on the parameters 
appearing in the Hamiltonian shared by the SACM and 
PACM, including also the x and y variables controlling 
the transitions between the three dynamical symmetries. 
The phase space was reparameterized in terms of three 
parameters A, B and C. 

In the case of the PACM, the potential energy sur- 
face was a sextic oscillator in the intercluster distance 
variable, while for the SACM the potential shape was 
more complex due to the restrictions enforced by the 
Pauli principle. The potential energy surface typically 
contained up to two minima, one spherical and one de- 
formed. The analysis identified both first- and second- 
order phase transitions for the PACM and the SACM, 
while in the latter case a critical line was also found. 

The results were illustrated with numerical studies on 
the 16 0+a and 20 Ne+a systems, which correspond to 
two spherical clusters and to one spherical and one de- 
formed cluster, respectively. The 5t/(3) limit was found 
to be the most appropriate one in reproducing the data of 
the cluster systems. Clear phase transitions were iden- 
tified in the x parameter controlling the transition be- 
tween the 5C/(3) and 50(4) limits. It was found that 
the PACM led to energy spectra that are rather different 
from the observed physical ones. 



Acknowledgments 

We gratefully acknowledge financial help from 
DGAPA, from the National Research Council of Mexico 
(CONACyT), OTKA (grant No. K72357), and from the 
MTA-CONACyT joint project. POH acknowledges very 
useful discussions with Octavio Castanos (ICN-UNAM), 
related to the differences in phase transitions in finite 
systems to the use of coherent states. The authors are 
also thankful to Jozsef Cseh for illuminating discussions 
on the subject. 



[1] H. Yepez-Martfnez, P. R. Fraser, P. O. Hess, and G. 

Levai, Paper I. 
[2] J. Cseh, Phys. Lett. B 281 (1992), 173. 
[3] J. Cseh and G. Levai, Ann. Phys. (N.Y.) 230 (1994), 165. 
[4] F. Iachello, Phys. Rev. C 23 (1981), 2778. 
[5] K. Wildermuth and Y. C. Tang, A Unified Theory of 

the Nucleus (Friedr. Vieweg & Sohn Verlagsgesselschaft 

mbH, Braunschweig, 1977). 
[6] P. O. Hess, G. Levai and J. Cseh, Phys. Rev. C54 (1996), 

2345. 



[7] O. S. Roosmalen, Algebraic Description of Nuclear and 

Molecular Rotation-Vibration Spectra, (Ph.D. Thesis, 

Groningen, Netherland, 1982). 
[8] O. S. Roosmalen and A. E. L. Dieperink, Ann. Phys. 

(N.Y.) 139 (1982), 198. 
[9] A. Leviatan and M. W. Kirson, Ann. Phys. (N.Y.) 188 

(1988), 142. 

[10] W. Greiner, L. Neise and H. Stocker, Thermodynamcis 
and Statistical Mechanics, (Springer, Heidelberg, 1995). 
[11] E. Lopez-Moreno and O. Castanos, Phys. Rev. C 54 



12 



(1996), 2374. 

[12] O. Castanos, R. Lopez-Pena, J. G. Hirsch and E. Lopez- 
Moreno, Phys. Rev. B 74 (2006), 104118. 

[13] R. Gilmore, Catastrophe Theory for Scientists and Engi- 
neers (Wiley, New York, 1981). 

[14] H. Yepez-Martmez, J. Cseh and P. O. Hess, Phys. Rev. 
C 74 (2006), 024319. 

[15] L. Parra Rodriguez, Master Thesis, UNAM, 2011. 

[16] P. R. Fraser, H. Yepez-Martmez, P. O. Hess and L. Parra- 
Rodriguez, J. Phys. Con]. Proc. (2011), accepted. 



[17] Cejnar P and Iachello F, 2007 J. Phys. A 40 581 

[18] Y. Zhang, Z.-F. Hou, H. Chen, H. Wei and Y.-X. Liu, 

Phys. Rev. C 78 (2008), 024314. 
[19] A. Frank, P. Van Isacker, Symmetry Methods in 

Molecules and Nuclei, (SyG editores, Mexico D.F., 2005). 
[20] G. Rosensteel and J. P. Draayer, Nucl. Phys. A 436 

(1985), 445. 

[21] R. Bijker and F. Iachello, Ann. Phys. 298 (2002), 334. 



